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Abstract. We report on the discovery of an unresolved X-ray 
. source inside the supernova remnant G3 15.4-2.3 (RCW 86). 

■ The source is located 7' to the Southwest of the geometrical 

■ centre and may be close to the actual explosion centre of the 
' supernova, which makes this a candidate for the stellar remnant 

associated with RCW 86. However, the presence of a possible 
optical counterpart with ^ 14 at 3" from the X-ray position 

. and evidence for long term variability means that the source 

' is probably an active star. A better X-ray position and better 
X-ray spectroscopy along with an identification of the optical 

' source are needed to exclude the X-ray source as a neutron star 

. candidate. 
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Table 1. A summary of ROSAT and Einstein observations 
which cover the point source. Dates refer to the start date of 
the observation. 



Date ID Target Instrument Exposure 

d/m/y s 



4/2/1980 


hl437s62 


SW 


Einstein HRI 


8756 


4/9/1992 


rp500078nOO 


NE 


ROSAT PSPC 


3212 


2/2/1993 


rp500079n00 


SW 


ROSAT PSPC 


4567 


3/2/1993 


rp500078aOI 


NE 


ROSAT PSPC 


2224 


18/8/1993 


rp500078a02 


NE 


ROSAT PSPC 


4509 


2/9/1994 


rh500077aOI 


SW 


ROSAT HRI 


2447 


4/8/1995 


rh500077a02 


SW 


ROSAT HRI 


6034 



1. Introduction 

In recent years it has become clear that young neutron stars do 
not necessarily manifest themselves as radio pulsars. Instead 
a large variety of unresolved objects associated with supernova 
remnants are thought to be the stellar remnants of the explosion 
(see Helfand 1998| for a review). Examples are the point-like 
sources recently discovered in Cas A (Tananbaum 1999 ) and 
Puppis A (Petre et al. 1996|), the enigmatic variable source in 



RCW 103 (Gotthelf et al. 1999), and a handful of relatively 



slow rotating X-ray pulsars called "anomalous X-ray pulsars" 
or AXPs (see Mereghetti 1998 for a review). 

Here we report on our analysis of an unresolved X-ray 
source in the supernova remnant RCW 86 (G3 15.4-2.3, MSH 
14-65). We discovered it duri ng our work on the X- ray pr oper- 
ties of the remnant (Vink et al. 1997 Bocchino et al. 200C ). The 



source qualifies as the possible stellar remnant associated with 
RCW 86, but the presence of a possible optical counterpart and 
long term source variability make this identification, as we will 
show, uncertain. 

RCW 86 is the ca ndidat e remn ant of the supernova AD 185 
(Clark & Stephenson |1977[ Strom |l99§), but the interpretation 
of the Chinese records is ambiguous (Chin & Huang 1994). 



The large extent of the remnant (40') can only be reconciled 
with an explosion as recent as AD 185, if its dista nce do es not 
exceed 1 kpc too much. However, Rosado et al. (1996) have 
pointed out that the kinematic distance towards RCW 86 seems 
to be higher, namely 2.8 ± 0.4 kpc. In that case RCW 86 may 
be physically associated with an OB association (Westerlund 



1969). The presence of a stellar remnant in RCW 86 would 
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establish the nature of the supernova as a core collapse super- 
nova (Type II or Ib/c), and therefore a connection with the OB 
association would be more likely. 



2. Data analysis 

Our analysis is based on archival ROSAT PSPC and HRI data 
and on an observation of RCW86 by the Einstein HRI. Both 
HRI instruments have a similar spatial resolution of 4" FWHM, 
but the ROSAT HRI is more sensitive. The PSPC has a resolu- 
tion of only 30" FWHM, is more sensitive than the HRI and 
has better spectral capabilities covering the energy range of 0. 1 
to 2.4 keV with a spectral resolution l\E/E = 0.4 at 1 keV. 
The data used for our analysis are summarised in Table |l| The 
starting point of the analysis are the basic screened event lists. 
Further processing, like photon extraction and barycentric cor- 
rection, was done with NASA's/fooZi v4.2 package. 
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Fig. 1. RCW 86 as observed by the PSPC (SW pointing). The 
circle indicates the position of the point source. 

Table 2. Best fit positions and count rates of the point source. 
The first three rows Hst the result of maximum likelihood fits 
to the ROSAT and Einstein HRI data. The statistic -2 In A has 
a distribution with three degrees of freedom. Position errors 
do not include systematic errors and correspond A(— 2 In A) = 
4.6 (90% confidence regions). The PSPC positions and count 
rates were estimated with a wavelet analysis method (Damiani 
et al. |1997[) using the energy channels 20-200 0.2 - 2 keV). 



ID 


RA (J2000) 


DEC (JOOO) 


error 


count rate 


-2 In A 




h m s 


o / // 




10~^cnts/s 




hl437s62 


1441 51.52 


-62 36 13.1 


2.6 


1.2± 0.4 


28.6 


rh500077a01 


14 41 51.64 


-62 36 11.9 


3.5 


2.1 ± 1.1 


12.4 


rh500077a02 


1441 51.11 


-62 36 13.2 


1.9 


3.4 ± 0.8 


55.3 


i-p500078n00 


14 41 53.3 


-62 36 37 


65 


10.9 ± 3.9 




ip500079n00 


1441 51.4 


-62 36 8 


15 


17.9 ± 2.7 




ip500078a01 


14 41 49.0 


-62 35 52 


59 


19.4 ± 6.5 




ip5{)0078a02 


14 41 48.9 


-62 36 22 


39 


26.8 ± 6.8 





2.1. Position of the source 

In order to find an accurate position for the unresolved source 
we used the HRI data and fitted the point spread function to the 
source using the very sensitive maximum likelihood fits (Cash 



1979 and e.g. Hasinger et al. 1994 1. We used a field of 80" x 



80" roughly centred on the source. The normalisation was cou- 
pled to the background level using iVtotai = ^bgd + ^src, 
-^totai being the observed number of photons and iVbgd, ^src 
the fitted background and source normalisations. The results of 
these maximum likelihood fits are listed in Table ^ which also 
lists the source parameters for the PSPC data as found with a 



wavelet analysis (Damiani et al. 1997). One of the HRI detec- 
tions has a significance level of 2.7 cr, not good enough to claim 
a detection by itself, but combined with the other detections the 
source parameters can be regarded as meaningful. 



Fig. 2. A 50" X 50" field image taken from the Digital Sky 
Survey 2. Overplotted is a circle with a radius of 5" centered 
on the point source and roughly corresponding to a 95% confi- 
dence region. 



Unfortunately, positions based on ROSAT sometimes suf- 
fer from errors in the attitude calculations which are typi- 
cally 6" (Hasinger et al. 



1994 and private communication). 



The observed scatter in the positions based on the Einstein and 
ROSAT HRI images, as compared to the statistical position er- 
rors, suggests that also our results are affected by systematic er- 
rors. Einstein HRI observations are less affected by systematic 
position errors, the typ ical sy stematic position error being ~ 2" 
(Van Speybroecket al. 1979| ). Adding statistical and systematic 
errors in quadrature the weighted average of the positions based 
on HRI images is a =14h 41m 51.42s and 5 = -62°36'12.9" 
(J2000) with a 1(7 position error of approximately 3". For a 
two dimensional gaussian this translates into a 95% confidence 
radius of 5". Note that a very bright, unresolved, radio source 
inside the remnant with coordinates a = 14''41™44.5'' and 
(5 = -62°34'47" (J2000) is clearly not as sociated with the un- 
resolved X-ray source (Dickel et al. 2000). 



2.2. Spectral analysis 

For the spectral analysis of the PSPC data of the unresolved 
source we extracted photons using a circular area with a ra- 
dius of 32" for the SW (on-axis) pointing and 44" for the other 
pointings. We estimate that with such radii we cover roughly 
90% of the photons coming from the unresolved source (c.f. 



Hasinger et al. 1992). Background spectra, extracted from an 
annulus around the source, were appropriately scaled and sub- 
tracted from the source spectra. The combined spectrum con- 
sists of 177 net source counts. The spectrum was rebinned in 
order to have at least 15 counts per bin. 



Jacco Vink et al.: An unresolved X-ray source inside the supernova remnant RCW 86 



3 



Table 3. Results of the spectral fits to the PSPC data. The luminosities are normalised to a distance of 1 kpc. For each model the 
spectrum was fitted with iVn — 1.7 10^^ cm^^ and with A^h as a free parameter Error ranges correspond to Ax^ = 2.7, or 
90% confidence limits. 



Model 


Parameter 


Normalization 






10^1 cm-2 


Lx(0.2-2.5keV) 
10=*^ erg/s 


X^/d.o.f. 


Blackbody 


kT 


= 0.17 ± 0.04 keV 


(4.6+JV) cm^ 






1.7 


3.6 


8.4/9 




kT 


= 0.23 ± 0.06 keV 


(0.59+J:L) cm^ 








1.6 


1.1/8 


Power Law 


r = 


o Q+0.12 
"-■••^-0.9 


(9.7 ± 2.7) lO^'^ ph/s/keV 


@ 


1 keV 


1.7 


16.9 


4.1/9 




r = 


2.9 ± 1.0 


(7.5 ± 2.9) 10^® ph/s/keV 


@ 


1 keV 


0.84 ± 0.47 


5.1 


1.7/8 


hot thin plasma 


kT 


= 0.68 ± 0.41 keV 


(6.8 ± 2.4) lO'^^ cm"^ 






1.7 


1.9 


14.7/9 


(mekal) 


kT 


= l-12tS:57 keV 








< 0.3 


1.3 


4.3/8 


idem,Z = 0.1 


kT 


= 0.63 ± 0.35 keV 


(4.4+*:^) 10^^ cm-3 






1.7 


3.0 


10.5/9 




kT 


= 0.87tUl keV 


(2.6 ± 0.7) lO'^" cm"^ 






0.3 ±0.2 


1.9 


1.0/8 



For our spectral analysis we us ed he spectral fitting pro- 
gram SPEXQ (Kaastra et al. |1996| ). Since we want to know 
whether the source qualifies as the potential stellar remnant 
associated with RCW 86, we fitted the spectrum with several 
emission models both with the interstellar absorption value 
fixed at Nn = 1.7 lO^'^^ cm^^, the typical absorption value 
for the X-ray emission of the supernova remnant (Vink et al. 



1997), and with A^h as an additional free parameter. The re- 
sults are listed in Table |[ The best fit values of A'h for all 
models seem to be in favor of a low absorption column towards 
the source, but also models with fixed A^h give acceptable re- 
duced values. The fact that models with three parameters 
result in very low reduced values (i.e. far from the expec- 
tation value), suggests that the statistics of the data is not really 
good enough to fit models with three or more parameters. All 
models provide reasonable fits to the data with only the thin 
plasma model with solar abundances and fixed iVn having a re- 
duced substantially larger than 1 . The spectrum appears to 
be rather soft as indicated by the steep power law index, F, and 
the low black body temperature. 



As for the variability on the timescales of month, at first 
sight there is little evidence for variability as all measured 
PSPC count rates are consistent with a count rate of (17.0 ± 
2.8) 10^'^ cnts/s. However, if we convert the Einstein and 
ROSAT HRI count rates to PSPC count rates using the best 
fit power law model in Table |^ (the conversion factors are 4.9 
and 2.7, respectively) we get the following PSPC count rates 
(in the same order as in Table ||): (5.8 ± 1.9) 10"'^ cnts/s, 
(5.7±2.7) 10-3 cnts/s, and (9.0±2.1) 10"^ cnts/s. The depen- 
dence of the ROSAT HRI/PSPC conversion factor on the cho- 
sen model is small (9%), but the conversion factor for the Ein- 
stein HRI count rates is more model dependent, varying from 
4.6 to 8.0. Even taking into account the model uncertainties it 
is clear that the observations are not consistent with a constant 
source count rate, although it is a strange coincidence that low 
count rates were only observed by the HRI instruments. Source 
contamination with the PSPC instrument seems unlikely, as no 
other unresolved sources are seen with the HRI instruments 
near the point source. Therefore, the X-ray source is very likely 
variable on a time scale of months to years. 



2.3. Timing analysis 

We searched the four PSPC observations for possible pulsa- 
tions using the Rayleigh method (Buccheri et al. 1983). This 



method is one of the most sensitive methods and it does not 
involve any binning of the data. A sensitive method is needed 
as the longest PSPC observations yielded only 101 events. We 
searched in each set for pulsations in the period range 0.02 to 
300 s, sampling the frequency range with step of 1 /Tobs with 
Tobs the total length of the observation. We compared the pe- 
riodograms to look for peaks showing up in two or more peri- 
odograms at or near the same period. Such correlations were, 
however, not found. The peak values of the Rayleigh statistic, 
Zf ^ 28, imply an upper limit to the pulsed fraction of ~20%. 



The black body model in SPEX vl.lO contains a small bug which 
we fixed for this analysis. The thin thermal plasma or CIE (collisional 
ionization equilibration) model is similar to the mekal model in xspec. 



3. Discussion 

So, could this X-ray source be the stellar remnant associated 
with RCW 86? The radius of the star, as inferred from the black 
body fit, is 1.7 km at a distance of 2.8 kpc. This is too small for 
a neutron star, but the spectrum may not be a black body. The 
X-ray luminosity is lower than the luminosity of typical AXPs, 
but it is consistent with the surface luminosity of young neutron 
stars in case pion cooling is important (Umeda et al. 2000). 



Also emission from a black hole accreting supernova fall back 
material should be considered as an alternative model (Umeda 



et al. 2000, Chakrabarty et al. 2000). The softness of the X- 



ray spectrum is a property shared with the unresolved sources 



in Puppis A, G296.5H-10.0 and AXPs (Mereghetti et al. 1996) 



Its long term variability is more in line with the behaviour of 
an active star, but a well established neutron star candidate like 
the point source in R CW 1 03 is also variable on similar time 
scales (Gotthelf et al. |l999|). 
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The main argument against an unambiguous identification 
of the source with a stellar remnant is the presence of a can- 
didate optical counterpart. The USNO-A2.0 catalogue (Monet 



0.5 1 1.5 2 

Energy (keV) 

Fig. 3. The PSPC spectrum of the unresolved source in RCW 
86. It consists of a combination of 4 individual observations. 
The solid line shows the best fit black body model with Nn 
fixed to a value of 1.7 10^^ cm^^. 



The position of the X-ray source is roughly at an an- 
gular distance of 7' from the geometrical center of the su- 
pernova remnant. This corresponds to a transverse velocity 

1200 km/s, if the supernova remnant is at a distance of 
2.5 kpc and 5000 yr old, or if the distance is 1 kpc and an 
age of 1800 yr. Such a kick velocity is rather high, but still 
consistent with observations of other neutron stars (Lyne & 
Lorimer 1994). Note, however, that the remnant is far from 



circularly symmetric, and the contrast in emission between 
the Northeast and Southwest of the remnant (see e.g. Vink 



et al. 1997, Bocchino et al. 2000l suggests that in the South 



west the shock wave is encountering a denser medium, which 
could mean that the actual explosion center was more to the 
Southwest of the geometrical center and closer to the unre- 
solved source. Interestingly, the remnant is quite symmetric 
in the Northeast/Southwest direction and the X-ray source is 
roughly located on the axis of symmetry. The position of the 
point source is one of its salient properties. The wavelet de- 
tection code has also detected another significant point source 
inside the remnant, 7' from the Northern shell with coordinates 
a ==14h 43m 47.0s and S = -62°19'29". The source seems em- 
bedded in a small 1 ' size extended structure. However, spectral 
analysis of this source rules out models with non-thermal or 
black-body radiation and the spectrum is consistent with ther- 
mal X-ray emission, characteristic for the North rim as reported 
by Bocchino et al. ( 200C| ). This suggests that this a small frag- 
ment of the shell seen in projection. 

It seems unlikely that the unresolved X-ray source is much 
further away than RCW 86, as the X-ray absorption column is 
too low for a distant source near the galactic plane. This means 
that it is unlikely that it is either an AGN or an X-ray binary. 
For the source to be a typical X-ray binary with a luminosity of 
Iq36 _ ]^q38 erg/s, it should be at a distance of 300 kpc. This is 
incompatible with both the X-ray absorption and with the size 
of our galaxy. 



et al. 1998) contains one source at an angular distance of 3" 
to the unresolved X-ray source (see Fig. The star has mag- 
nitudes r = 13.6 and b = 15.9, which suggest a K-star, al- 
though the magnitude errors and the unknown extinction make 
a late type G-star or a early type M star also possible (Zombeck 
1990). The typical V magnitude of this star would heV^ 14, 



which, for a main sequence star, implies a distance of 250 pc. 
For this distance the X-ray luminosity of the X-ray source is 
Lx ^ 6 10^^. This is consistent with the X -ray luminosity of 
a bright active star (see e.g. Zombeck 1990 p. 88, Agrawal et 
al. 1986). The possibility that this is just a chance alignment 



can be estimated using the USNO-A2.0 catalogue. There are 
about 0.0015 stars per square arcsec with ttir < 16 in the field 
around the X-ray source, which means that the chance to find a 
bright star within the 95% error circle (see Fig. ||) is ^11%. 

The evidence for long term variability and the presence of 
a late type star close to the X-ray position favour the identifica- 
tion of the X-ray source with an active star, rather than with the 
stellar remnant associated with RCW 86. However, to firmly 
establish this identification, optical spectroscopy of the pos- 
sible optical counterpart and X-ray observations with a better 
positional accuracy of the source are needed. Improved X-ray 
spectroscopic data can also be used to distinguish between the 
typical optically thin thermal spectrum of an active star and the 
featureless spectrum of a neutron star. 
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